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Abstract 

We summarize our renormalization group approach for the vector 
model as well as the matrix model which are the discretized quantum 
gravity in one- and two-dimensional spacetime. A difference equation is 
obtained which relates free energies for neighboring values of N . The 
reparametrization freedom in field space is formulated by means of the 
loop equation. The reparametrization identities reduce the flow in the 
infinite dimensional coupling constant space to that in finite dimensions. 
The matrix model gives a nonlinear differential equation as an effective 
renormalization group equation. The fixed point and the susceptibility 
exponents can be determined even for the matrix models in spite of the 
nonlinearity. They agree with the exact result. 

1 Introduction 

The long-standing challenge for the quantum theory of gravity has met at least a partial success 
in low dimensions. The two-dimensional quantum gravity is particularly interesting from the 
viewpoint of string theories too. The discretized approach such as matrix models offers a pos- 
sibility for a nonperturbative treatment Exact solutions of the matrix model have been 
obtained for two-dimensional quantum gravity coupled to minimal conformal matter with central 
charge c < 1. It has been very difficult to obtain results for two-dimensional quantum gravity 
coupled to conformal matter with central charge c > 1. Although one can easily write down 
matrix models for cases with c > 1, these models are not solvable up to now^. The numerical 
simulations suggest that it is not at all obvious if a matrix model candidate to describe a c > 1 
model has the continuum description 0. Therefore it is useful to obtain approximation schemes 
which enable us to calculate critical coupling constants and critical exponents for unsolved matrix 
models, especially for c > 1. In order to make use of such a scheme, we first need to make sure 
that the approximation method gives correct results for the exactly solved cases. 

"The speaker. 



Brezin and Zinn- Justin have proposed a renormalization group approach to the matrix model 
Consequences of their approach have been examined by several groups 0. A similar approach 
has been advocated previously for the expansion in somewhat different contexts P]. The 
vector model has been proposed for a discretized one- dimensional quantum gravity, in the same 
way as the matrix model for a discretized two-dimensional quantum gravity [^.Recently we have 
analyzed the vector model by means of the renormalization group approach and have clarified 
its validity and meaning [|10].More recently, we have succeeded in extending our analysis of the 



renormalization group approach to matrix models [jTT . 

The purpose of this paper is to present a new simplified derivation of our results on vector 
models and to summarize some of our results on matrix models in ref. [|lT|. For the vector 
model, a difference equation is obtained relating the free energy — log Z]\f_2{g) to the free energy 
— log Z]\f{g — 26g) with the coupling constant shifts 6gk of order 1/N in infinitely many coupling 
constants. We obtain infinitely many identities which express the freedom to reparametrize 
the field space. Thanks to these identities, we can rewrite the flow in the infinite dimensional 
coupling constant space as an effective flow in the space of finite number of coupling constants. 
In this note, we shall simplify the derivation of the reparametrization identities by using a new 
method similar to the loop equation for the matrix model. The resulting effective beta function 
determines the fixed points and the susceptibility exponents which agree with the exact results 

We also find that a similar procedure works for matrix models. Namely we can obtain a 
difference equation relating the free energy — log Zjv_|_i(5f) to the free energy — log Z^^g + Sg) with 
shifts 5g of order 1/iV in infinitely many coupling constants. We find infinitely many identities 
expressing the reparametrization freedom of matrix variables. The identities are sometimes called 
loop equations 0, [|1^ and enable us to rewrite the flow as an effective renormalization group flow 
in the space of finite number of coupling constants. A new characteristic feature of the matrix 
model is that the resulting effective renormalization group equations is nonlinear contrary to 
the vector model. In spite of the nonlinearity, we can obtain the fixed points and the critical 
exponents |11|. As an explicit example, we analyzed the one matrix model with a single coupling 



constant and find a complete agreement with the exact results for the first critical point m = 2 

ra. 



2 Renormalization Group Approach for Vector Models 

It has been proposed that the free energy F{N, g) of the matrix model should satisfy the following 
renormalization group equation [B] 



F(N,a) = r(a), 



where I3{g) is called the beta function. The anomalous dimension and the inhomogeneous term 
are denoted as 7(5') and r{g) respectively. A fixed point g^, is given by a zero of the beta function. 
In the following we shall show that this renormalization group equation is valid for vector models. 

The partition function of the 0{N) symmetric vector model is given by 
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where is an iV dimensional real vector |p . Here we introduce infinitely many coupling constants 
Qk, since we need all possible induced interactions after a renormalization group transformation 
even if we start with a few coupling constants only. The expansion of the logarithm of the 
partition function gives contributions from h loops as terms with N^~^. The vector model has 
the double scaling limit N ^ oo with N^/^'^[g — g^) fixed, where the singular part of the free 
energy satisfies the scaling law P] 
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By inserting this form, one find that the critical point corresponds to the fixed point l3{g^,) = 
and that susceptibility exponents 70,71 are related to the derivative of the beta function 



71 = 7o = 2-7((7,)//?'((7,). 



(4) 



In order to obtain a difference equation, we start with the partition function Z7v_2(fl')- After 
integrating over angular coordinates in R^~^, we perform a partial integration in the radial 
coordinate x = d)^ 



ZN-2{g) 



^-1 

TT 2 

rxf^i) Jo 

TT 2 



exp 



9k k 



-(iv-2)Ei 



k=l 



X 



iN-2)Y, 



k=l 



2k 



(5) 



Identifying the right hand side with ((A^ — 2)gi/27T)ZN{g — 26g), we obtain a difference equation 
for the logarithm of the partition function 
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= -[{-\ogZ^{g-26g))-{-\ogZr,{g))]. 
where the shifts 6gk of the coupling constants are found to be 
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We stress that no approximation is employed to obtain this difference eq. apart from neglect- 
ing the higher order terms in contrary to the perturbative method advocated in ref.[^]. 

One can bring the quadratic term in the potential to the standard form 0^/2 since gi can be 
absorbed by a rescahng gicj)'^ — > 0^. We have 

^Tvlfi-i, 5^2, 5-3, ■••)= 5'i~^^^^7v(l, ^2,^3, •• •)> 9k = 9k/gi- (8) 
We shall define the free energy for the vector model 



F{N, gi,g2,g3, 
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We sometimes write the free energy as a function of g2, gs, ■ ■ ■, since it is independent of gi because 
of the rescahng identity (|) if one uses gi,g2,---as independent variables. We denote the partial 
derivatives with respect to gi, ^2, ^3, • " " by \g, and those with gi, g2, gs, - ■ ■ by \g 
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In the N 00 limit, we can obtain a differential equation from the difference equation 
as a renormalization group equation for the free energy F 
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This renormalization group equation clearly shows a flow in the infinite dimensional coupling 
constant space as is usual in the Wilson's renormalization group approach. Eq.(|TT]) shows that 
the anomalous dimension is given by 7(^) = 1. Therefore two susceptibility exponents are related 
by 7o + 7i = 2. 



3 Reparametrization Identities 

Our key observation in using the renormalization group equation is the ambiguityto identify the 
renormalization group flow in the coupling constant space [0]. Though the above equation (^) 
seems to describe a renormalization group flow in the infinite dimensional coupling constant 
space, the direction of the flow is in fact ambiguous because all the differential operators {d/dgk) 
are not linearly independent. Since the model is 0{N) invariant, we can obtain new informations 
only from reparametrizations of the radial coordinate x = (f)^ which become reparametrizations 
of a half real line. In order to facilitate the derivation of the reparametrization identities, it is 
useful to make the following reparametrization 
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(12) 



where e is an infinitesimal parameter and C is a parameter to generate all possible reparametriza- 
tions as a power series. Substituting (|T2p in (|^), we obtain an identity 

= 0, (13) 
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we shall call this identity a loop equation since it resembles that in the matrix model |^. It is 
convenient to define the expectation value W{C) of a resolvent 



C-ct>^ C^+i C ±^.0+^ do, 
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By expanding the identity (|T3D, we can obtain various reparametrization identities which form a 



representation of the Virasoro algebra |]T0|. We see that derivatives of the free energy in terms of 
infinitely many coupling constants Qk are related by infinitely many reparametrization identities. 
Thus one can expect that only a finite number of derivatives are linearly independent. The loop 
equation for the vector model (|T3p can be solved as 



1 - 2CV'{0 
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(16) 



The resolvent to the leading order in expansion is given by omitting the last term. 

To illustrate the use of the reparametrization identities, we shall first take the case of a single 
coupling constant. Let us consider a point in the coupling constant space {91,92,93^94:: ■■ ■) = 
{91, 92, 0, 0, . . .). The coupling constant shifts 69k given by eq.(0) are then found to be 
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Inserting the shifts 59k to the right hand side G of the renormalization group equation ([TTD, we 
find 
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By using the solution for the resolvent (|T6D to the leading order in 1/A^, we find that the renor- 
malization group flow in infinitely many coupling constant space is now reduced to an effective 
flow in single coupling constant 92- Thus we obtain a renormalization group equation with the 

np 

effective beta function P (92) and the inhomogeneous term r(g2) 



ON 



092 



(19) 



where the effective beta function and the inhomogeneous term are given in terms of A = 4^2 + 1 
as 
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The effective beta function /3^^(^2) exhibits a zero at ^2 = = 92*- Furthermore, we get 

the susceptibility exponents 71 = 3/2, 70 = 1/2 from the derivative of the effective beta function 
by using @) . The fixed point and the susceptibility exponents are in complete agreement with 
the exact results for the m = 2 critical point of the vector model corresponding to pure gravity [0]. 
The beta function P {92) has also a trivial fixed point at 92 = 0, which is ultraviolet unstable 
since d(3^^/dg2{g2 = 0) = -1 < 0. 



By a systematic expansion of the free energy in powers of 1/iV, we can extract the complete 
information from the renormahzation group flow namely from the difference equation and the 
reparametrization identities 

oo 

F(iV,^2) = E^''A(^2). (22) 

h=0 

We should evaluate the difference equation (|]) and solve the loop equation ( |16|) both up to the 
order N~^. Then we obtain an ordinary differential equation for fh 



(1 - h)M~g,) - /3^^(~g,)^(g,) = r,{~g,), (23) 

dg2 

where the effective beta function is common to all h, but the inhomogeneous terms depend 
on h. We see immediately that the general solution is given by a sum of an arbitrary multiple 
of the solution of the homogeneous equation and a particular solution of the inhomogeneous 

np 

equation. Since both the effective beta function /3 (^2) and the inhomogeneous term r(^2) are 
analytic in ^2 around the fixed point ^2*, the singular behavior of fh comes from the solution 
of the homogeneous equation. It is important to notice that the singular term corresponding to 
the continuum physics (the so-called universal term) is specified by the beta function alone. The 
normalization ah of singular terms, however, cannot be obtained from the renormahzation group 
equation 

fhih) = fh{,g2)sing + fhih) analytic 
U92)sing = ih - ^2*)^'-^^)"a, + ■ ■ ■ . (24) 

If we sum over the contributions from various h, we find that the renormahzation group equation 
determines the combinations of variables appropriate to define the double scaling limit. On the 
other hand, the functional form of the singular part of the free energy Fging in the scaling variable 
N^^^^{g2 — 92*) is undetermined corresponding to the undetermined normalization factor ah for 
the singular terms of each h 

oo oo 

F(iV, ~g)s^ng = E N''^ M~92) s^ng = E ^"'(^2 - ^2*)'~"°'"^'a, 
h=0 h=0 

= {92 - 92*?-'''f{N'/''^{92 - ^2*))_,. (25) 

Let us note that the free energy is completely determined if we impose the condition F{N,g2 = 
0) = 0, which follows from the definition (^. 

We can extend the analysis to more general situations of finitely many coupling constants, 
such as the case of the two coupling constants gi = 1, ^27^ 0, ^3 7^ 0, ^fc = (A; > 4). 
We can write down the shifts of coupling constants 6gk in eq.(|^) explicitly, and we can obtain 
the resolvent as a solution of the loop equation (0). By inserting these to the right hand side 
G of the renormahzation group equation, we obtain effective beta functions and inhomogeneous 
terms. As reported in ref. [0, we found a perfect agreement with the exact result for m = 2, 3 



cases. Moreover, we were able to explicitly draw the renormahzation group flow diagram in the 
space of two coupling constant space [|l^]. Higher multicritical points can be analyzed similarly. 



4 Nonlinear Renormalization Group Equation for Matrix 
Models 



In this last section we report briefly on our new results on matrix models |TT|. The partition 
function Z^lg) of the matrix model is defined by an integral over an N x N hermitian matrix $ 
with a potential V{^) 



ZN{g) = /d^'$exp[-iVtrV^($)], V{^) = Y.^-^^' 

k=l 



(26) 



The cubic interaction with a single coupling constant g corresponds to V^($) = + We 
can integrate over the angular variables to obtain an integral over the eigenvalues A M 
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where denotes the Van der Monde determinant AAr(A) = ni<i<j<Af(-^i ^ cind cat 
In order to relate Z]\r^i to Zjv, we shall integrate one of the eigenvalues \n+i in ^a^+i 
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The Aat+i integral can be evaluated by a saddle point method to the leading order in The 
saddle point A^ , ^ is determined as a function of tr by the saddle point equation 
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By inserting the saddle point X%^i into the partition function, we obtain 

Zj,+,{g) ^ c^+i /^^''^exp [-{N + 1) tr V($) - NVjXl,^,) + 2tr log [A^^^ - + Q(iV°) 
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Here ( ) represents the normalized average with respect to the measure exp [— A^tr V^($)]. 
We observe from the above equation and eq.(^) that infinitely many numbers of operators of 
the form tr $™ tr $" ■ ■ ■ are induced. 



The situation, however, simphfies if we appeal to the large- limit. In this limit a multi-point 
function of f7(iV) -invariant operators factorizes into a product of one-point functions 



{O 0'---) = {0){0')--- + 0[ — ]. 



(32) 



In favor of this factorization property, eq.(0) reads 
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where (A^, ^) is determined in terms of (tr^™) by averaging eq.(|30D 
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(34) 



Let us introduce the free energy which is normalized to vanish for vanishing coupling constants 
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where Zn^Qi = 0, g2 = I, gk = (k > 2)) = 2^/^{7c /N)^^/^ . By taking the ^ oo hmit, we find 
the following differential equation as a renormalization group equation for the free energy 
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We observe that this renormalization group equation describes a flow in infinite dimensional 
coupling constant space analogous to the vector model. We stress that this RG equation does 
not involve a perturbation with respect to coupling constants, unlike the approximation schemes 
proposed in ref. 0. 

Similarly to the vector model, we can find infinitely many identities expressing the reparametriza- 
tion freedom of matrix variables. In order to derive these identities, it is most convenient to 
formulate in terms of the loop equations 0], |12|. Let us define the expectation value of the 
resolvent 

"'■K' = {"(c^)) P"' 



By a reparametrization of matrix variables, we can obtain the loop equation to the leading order 
in 1/iV explicitly 



wiO'-v'iOWiC) + QiC;V) = o, 



(3J 



If we consider the space of single coupling constant for the cubic interaction {gi, g2, gs, ■ ' ') = 
(0, 1, (?3, 0, ■ ■ ■), we can relate derivatives with respect to all the other coupling constants to that 
with rpsnpp.t tn thfi sinprlp, rnimlinp" r.nnstarit rin bv mparis nf thp, Innn pniiatinn 



Therefore we obtain to the leading order in 

. ). m 

Contrary to the vector model, this effective renormalization group equation contains nonlinear 
terms in dF^/dgs. 

This nonlinearity is the most characteristic feature of the matrix model. 

In spite of the nonlinearity, we can find the fixed points and the critical exponents in the 



following way |]Tl|. If we expand the coefficients Pn in powers of the coupling constant around 
the fixed point (73*, they are analytic 

00 

Pn{g3) = T.(^nk{g3-93*)'- (41) 
k=0 

The non-analyticity of the free energy should come from solving the differential equation 
We assume that the free energy has singular and analytic terms as in the vector model (|2 

oo oo 

Fi^^ = Fanalytic + Fging, Fanalytic = X! ^^{93 — 93*)'' 1 Fging = ^ dk{93 — 93*)^^'^ ■ (42) 
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By comparing the power series expansion of the renormalization group equation, we find consis- 
tency conditions for the above expansion to be valid. The first four conditions read 
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By solving these four equations we can determine four quantities, namely the fixed point, the 
coefficient ci, the critical exponent 7, and the coefficient C2. Therefore these equations are enough 
to determine the fixed point and the critical exponents. We have found that the above set of 
equations in fact has a solution that agrees with the exact result [11| g^^ = ±2~^ -3"^/^, 71 = 
7 = 5/2. 

Similarly we can determine all the coefficients Ck {k > 0) and dk/do {k > 1) except for the 
overall normalization of the singular term do. Moreover we can show that higher genus free 
energies have singular terms precisely needed for the double scaling behavior 

00 

F{N, gs)sina = E A^"'V/.(^73).^n„ fh{93)s^na = (93 - 93*Y'-''^^'d'^ + ■■■■ (44) 



The singular terms are determined up to the normahzation for each genus. These features 
are analogous to the case of the vector model. However, a distinctive feature of the matrix 
model is that the two coefficients ci and C2 are needed to determine the fixed point and the 
critical exponent 71. This novel new feature is a direct consequence of the nonlinear nature of the 
renormalization group equation for the matrix model (0). It is intriguing to observe that these 
two coefficients can be determined self-consistently without referring to all the other coefficients 
Ck {k 7^ 1, 2). We can perform a similar analysis for cases with more couplings corresponding to 
the higher multicritical points, and also multi- matrix models [pJ]] . 
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